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p ^ Abstract 

The Booster Neutrino Experiment (MiniBooNE) searches for oscillations 
using the ~lGeV neutrino beam produced by the FNAL Booster synchrotron. 
CN The array of photomultiplier tubes (PMTs) lining the MiniBooNE detector records 

^ Cherenkov and scintillation photons from the charged particles produced in neu- 

p\J trino interactions. We describe a maximum likelihood fitting algorithm used to 

reconstruct the basic properties (position, direction, energy) of these particles from 
the charges and times measured by the PMTs. The likelihoods returned from fit- 
^ ting an event to different particle hypotheses are used to categorize it as a signal Ve 

Q\ event or as one of the background processes, in particular charged current quasi- 

elastic scattering and neutral current vr'^ production. The reconstruction and event 
>• selection techniques described here can be applied to current and future neutrino 

experiments using similar Cherenkov-based detection. 
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1 Introduction 



Cherenkov detectors of >1 kiloton have played a key role in the estab- 
hshment of the phenomenon of neutrino oscillations [la-b,2a-c,3a-c]. These 
detectors typically consist of a large volume of a homogeneous transparent 
medium (water or mineral oil) with a high index of refraction surrounded by 
an array of photomultiplier tubes (PMTs). Cherenkov photons produced by 
charged particles emerging from the neutrino interactions and traversing the 
medium with (5 > 1/n (where n is the index of refraction and (5 — v/c) are 
detected by the PMTs. The photons are emitted at an angle 9c relative to the 
track direction, where cos 6c = l/nf3. The radiation is azimuthally symmetric 
about the track direction, resulting in a ring-like pattern that can be identi- 
fied on the PMT array. The quantity, spatial distribution, and arrival times 
of these photons provide information on the location, direction, and energy of 
the particle. We refer to the extraction of such information from the charge 
and time measurements of the PMTs as "reconstruction." An analysis of the 
ring profile can also provide information on the identity of the particle and, if 
multiple rings are detected, the number of particles. 

In this paper, wc discuss the event reconstruction algorithms used in the 
Booster Neutrino Experiment (MiniBooNE), which searches for an excess of 
i/g interactions indicative of — > oscillations using the Fermilab ~1 GeV 
neutrino beam. The MiniBooNE detector is a sphere of radius 610.6 cm filled 
with Marcol 7 mineral oil {n 1.47) and divided into two optically isolated 
regions by an opaque shell of radius 575 cm, concentric with the sphere. The 
surface of the inner "main" region is instrumented with an array of 1280 
inward-facing 8 in. PMTs which detect the light produced by the neutrino 
interactions. The outer "veto" region is instrumented with 240 PMTs and is 
used to tag charged particles that enter the detector from outside {e.g., cosmic 
muons) or that exit the main region. Though no scintillator was added to the 
Marcol 7, it scintillates weakly, resulting in the production of delayed isotropic 
light for particles with sub- and super- Cherenkov velocities. A more detailed 
description of the experiment can be found in Refs. [4] and [5]. 



2 Neutrino Interactions at MiniBooNE 



The ~lGeV neutrino beam at MiniBooNE results in interactions with 
relatively low outgoing particle multiplicity. The largest interaction channel is 
the charged current quasi-elastic (CCQE) process: 

ue + n^r+p (1) 
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which accounts for ~40% of all neutrino interactions in the MiniBooNE de- 
tector. Since the recoil proton of Eq. (1) is typically below Cherenkov thresh- 
old, only the outgoing lepton produces significant light. In modeling such 
events within a reconstruction algorithm, one can therefore consider only the 
Cherenkov and scintillation light produced by the outgoing lepton. While the 
recoiling nuclcon can produce scintillation light, this additional source of light 
is not considered in the reconstruction. 

A muon in the MiniBooNE detector exhibits minimum-ionizing behav- 
ior through most of its path with little chance of radiative energy loss. In 
contrast, electrons typically induce electromagnetic showers, resulting in ad- 
ditional electrons and positrons that emit Cherenkov light. Thus, muons and 
electrons have significantly different Cherenkov ring patterns. This difference 
is the basis for distinguishing a muon track from an electron track, and the 
reconstruction algorithm has a model for each. 

The next most abundant process in MiniBooNE is single pion production, 
which occurs primarily via A resonance or coherent scattering: 

+ TV ^ £- + iV' + TT (CC tt) 

(2) 

N ^ N' + 7: (NC tt) 



where A^^'^ denotes a nucleon in the case of resonance production and a nucleus 
in the case of coherent production. Of particular interest is neutral current 
(NC) 7r° production. The two photons from 7r° — > 77 decay induce electro- 
magnetic showers indistinguishable in MiniBooNE from those induced by elec- 
trons. Since the recoil nucleon/nucleus is typically below Cherenkov threshold, 
these NC 7r° events are well-described in the reconstruction algorithm by two 
electron-like tracks pointing back to a common vertex. The presence of two 
distinct Cherenkov rings allows for the separation of tt" events from electron 
{ve CCQE) events. However, vr'^ misidcntification (for example, due to a large 
energy asymmetry between the two photons or due to a small photon open- 
ing angle which leaves the Cherenkov rings overlapping) accounts for most of 
the reducible background in the — > oscillation search. Thus, successful 
reconstruction and identification of 7r° events is critical. 

While other channels can be accommodated by the reconstruction algo- 
rithm, the z/g appearance oscillation analysis uses only four event models: single 
electron track, single muon track, two 7 tracks, and two 7 tracks with a 7r° 
invariant mass. 
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3 The Detector Response 



The Marcol 7 mineral oil used in MiniBooNE exhibits a rich array of opti- 
cal phenomena despite its ~20 m extinction length near the peak of PMT 
sensitivity (~400 nm). Cherenkov and scintillation light production is ac- 
companied by photon absorption, fluorescence (with several possible excita- 
tion/emission spectra and lifetimes), and Rayleigh and Raman scattering. The 
left plot in Figure 1 summarizes the rates of various processes as a function of 
wavelength. The cumulative extinction rate is shown as the black line. In the 
near-ultraviolet region (<320 nm), fluorescence processes dominate the extinc- 
tion rate. The fluorescence lifetimes range from 1 — 35 ns. In the visible region 
(>320nm), the dominant processes are Rayleigh scattering and absorption. 




Fig. 1. Left: Rates of optical processes in Marcol 7 as a function of wavelength. 
The solid black line is the overall extinction rate obtained from laboratory measure- 
ments. The dashed black line is the extrapolated rate based on in situ data. The 
curves labeled "Fluor n" are the excitation rates for the four identified fluorescence 
processes. The light blue points represent the measured rates of Rayleigh scatter- 
ing, and the dashed light blue and gray lines represent extrapolated rates. Right: 
Reconstructed photon arrival times for R1408 (top) and R5912 (bottom) PMTs for 
397 nm light flashed from the center of the detector. The black histogram is the 
distribution from data and the blue is the complete Monte Carlo simulation. The 
green (red) shows the simulation with reflections (and scattering) turned off. 



The 1520 8 in. PMTs in MiniBooNE are of two types: 322 model R5912 
PMTs and 1198 model R1408 PMTs, both from Hamamatsu. AU of the R5912 
PMTs are located in the main PMT array. The PMT time response (particu- 
larly the late-pulsing behavior) and the variation of PMT efficiency with inci- 
dent angle have been characterized in external measurements using a pulsed 
LED [6] . The PMTs in situ have been studied using a laser calibration system 
that produces sub-nanosecond bursts of nearly isotropic light in the detector. 
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The right panel of Figure 1 shows the times of PMT hits recorded in laser 
calibration events relative to their expected times. The colored histograms 
come from Monte Carlo simulation, with the red curve showing PMT and 
electronics timing effects (including pre- and late-pulsing) , the green curve in- 
cluding photon scattering, and the blue curve adding photon reflections from 
the PMTs and the surface of the main detector region. The resulting time 
structure matches well with that seen in data (black points). We note that 
only the earliest photoelectron's time is reported when multiple photoelec- 
trons are present in a hit. (A fresh hit can be seen after an electronics dead 
time of ~300 ns.) The reported charge, however, does account for additional 
photoelectrons. 

Absorption, scattering, reflections, and fluorescence processes influence 
both the topology and the time structure of hits. Additionally, scintillation 
light, emitted with a lifetime of ~35ns, adds intrinsically delayed hght to 
the prompt Cherenkov photons. The event models within the reconstruction 
algorithm account for all these phenomena. In the discussion, we refer to 
Cherenkov and scintillation light that travels undisturbed from the particle 
track to the PMT as "direct" light. Light that undergoes fluorescence, scat- 
tering, or reflection we term "indirect" light. 

A detailed GEANT3-based [7] Monte Carlo simulation of the MiniBooNE 
detector serves as the central tool for developing the reconstruction algorithm's 
predictive models. The accuracy of the models influences performance through 
resolution and biases in extracted track parameters and through the particle 
separation capabilities of the likelihood-based hypothesis testing. To balance 
the need for likelihoods that are as accurate as possible with the desire to 
contain the computational demands, detector and track model information 
extracted from the Monte Carlo simulation is tabulated, approximated, or 
parametrized wherever possible while minimizing any sacriflce in performance. 
Also, thanks to the excellent stability of the detector's optical properties, a 
single configuration of the reconstruction is sufficient for use on all MiniBooNE 
physics data. 



4 The Single Track Model 



The reconstruction of a single track in the detector, whether an electron 
or muon, is based on a model with seven parameters: 

• the starting point: XQ,yo,Zo 

• the starting time: to 

• the direction: 6*0, 4>o 

• the kinetic energy: Eq. 
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The starting point and time are referred to as the "vertex" of the event. The 
direction is described in terms of the polar angle 6q defined with respect to the 
direction of the neutrino beam and the azimuthal angle 0o about this beam 
axis. For the electron model, the stochastic variations in the electromagnetic 
showers are modeled in average - no attempt is made to account for event- 
by-event variations in the shower profile. This is also true of the much smaller 
variations in the muon propagation {e.g., straggling). This simplification allows 
the track to be fully specified from its initial properties summarized in the 
seven parameters above. We refer to this vector of parameters as x. 

The observables of an event are the measurements from the 1280 PMTs 
in the main region of the detector. For each PMT, we have: 

• a bit indicating whether the tube registered a hit 

• if the tube was hit, the measured charge of the hit 

• if the tube was hit, the measured time of the hit. 

Assuming that the PMTs behave independently, the likelihood for an observed 
set of PMT measurements given track parameters x can be expressed as: 

£(x) = n (1 - hit;x)) X n hit;x) /,(g,;x) /i(t,;x) (3) 

unhit hit 



where the products are taken over all unhit and hit PMTs and: 

• P{i hit;x) is the probability that the ith tube is hit, given parameters x. 

• qi,ti are the measured charge and time on the ith PMT. 

• fq{qi;x) is the probability distribution function (PDF) for the measured 
charge on the ith PMT, given x, evaluated at Qi. 

• ft{ti;x) is the PDF for the measured time given x, evaluated at ti. 

It is convenient to work with the negative logarithm of C, and since the charge- 
and time-related portions of the likelihood are distinct, it is helpful to define: 

F(x) = -log£(x) = F,(x) + Fi(x) (4) 



where 



F,(x) = - ^ log(l - P{t hit;x)) -^log(P(^ hit;x)/,(g,;x)) 

unhit hit 

F*(x) = -^log(/*(i,;x)) 

hit 



For brevity, we refer to Fq and Ft as the charge and time likelihoods, respec- 
tively, although they are the negative logarithms of the hkehhoods. 
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4-1 The Charge Likelihood 



If the number of photoelectrons n produced in a PMT is known, then fq{q) 
for the observed charge q is fully specified in terms of the response properties 
of the PMT (assumed fixed) without any reference to any other property of 
the detector. Further, n is a Poisson variable whose mean is a function of 
the track parameters x. These considerations motivate a two-step approach 
to calculating the charge likelihood. 

In the first step, using the known optical photon and particle propaga- 
tion properties of the detector, one determines for a given particle type (elec- 
tron/muon) and set of track parameters x the average number of photoelec- 
trons that a particular PMT should observe. This quantity is referred to as //j, 
the "predicted charge" for the ith PMT. In calculating the predicted charge, 
one considers the quantity and angular distribution of light produced by the 
track, the absorption, scattering, and fluorescence of light as it traverses the 
mineral oil, the acceptance of the PMT, and anything else that influences the 
mapping x i— > {/li}. 

For the second step, Fq given in Eq. (5) can be rewritten in terms of the 
predicted mean charges /if 

F,(x) = -J2 log(l - P{i hit; //O) - E log(^(^ hit; //,)/,(?,; (6) 

unhit hit 



where P and fq are now functions of the {/ij}, with the mapping x i— > {//j} from 
the first step left implicit. These two functions depend only on the properties 
of the PMT response and are, in particular, independent of the happenings 
within the detector volume. This stepwise approach decouples the track and 
optical photon model from the PMT response model. The latter is established 
once using laser calibration data. The laser system provides a source of ap- 
proximately isotropic and prompt light with controllable intensity. The value 
of /i for a particular laser intensity can be determined from the occupancy of 
the PMTs. The function fq{q] n) for each PMT type can be determined by 
obtaining the distribution of measured charges q seen at various intensities 
These distributions incorporate all processes associated with the charge 
response, such as the cascade of charge down the dynode chain, the amplifi- 
cation and digitization of the anode signal, and the conversion of the digital 
readouts into calibrated charge values. Examples of fq{q;n) obtained in this 
way are shown in Figure 2. 
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R1408 PMTs 



1 23456789 10 
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Fig. 2. Examples of f{q; ji) for ji = 0.7 photoelectrons. Note that the distribution 
only includes cases where the PMT registers a hit. 

4-2 Calculating the Predicted Charge /i 

We begin our description of the predicted charge calculation with a sim- 
ple scenario involving a fictitious point source of isotropic direct light. This is 
then generalized to a line source of isotropic light, appropriate for the scin- 
tillation production of a track. We then consider a line source in which the 
emission has a non-trivial (but azimuthally symmetric) angular distribution, 
representative of Cherenkov radiation. Finally, we account for the contribution 
of indirect light {i.e., scattering, reflections, etc.) from both the scintillation 
and the Cherenkov emission. In order to simplify notation and without loss of 
generality, the discussion involves fixed but arbitrary PMT location and track 
parameters x. 

One needs only two parameters to describe the geometric relationship 
between an isotropic point source and a PMT. We choose (1) the distance 
r between the source and the PMT and (2) the angle of incidence rj at the 
PMT, with 77 = corresponding to incidence parallel to the PMT axis. The 
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Fig. 3. Geometry of a line source of isotropic scintillation light (left) and directional 
Cherenkov light (right). In each case s is the distance from the origin of the track, 
r(s) is the distance to the PMT from the point along the track, and t/(,s) is the angle 
of PMT incidence for light emitted from this point. The presence of Cherenkov light 
requires the additional parameter 9{s) describing the light emission angle. 

predicted charge on the PMT can be written 

Aipoint.sci = $sci ^{r) Tsci(r) e{ri) (7) 



where <l>sci is an event-energy-dependent scintiUation hght yield {i.e., $sci = 
$sci(-^o) ); ^(^) is a r-dependent sohd angle factor, Tsci(r) is the transmission 
of the oil and PMT glass as a function of light propagation distance, and €{rj) 
is the acceptance of the PMT as a function of the angle of incidence. The 
transmission depends on the wavelength spectrum of the light source and, 
thus, is specific in this case to scintillation. The wavelength-dependent PMT 
photocathode efficiency is also included in T(r). The sohd angle factor and 
PMT angular acceptance are purely geometric, independent of the particulars 
of the light production. 

To generalize to an extended source of scintillation light, as depicted in 
Figure 3(left), the latter three factors in Eq. (7) become functions of s, the 
distance along the track from its origin. We must also include an emission 
profile Psci{s) describing the distribution of scintillation production along the 
track. This profile satisfies /^p(s)rfs = 1. The predicted charge can now be 
obtained by integrating along the track's axis: 

oo 

A«sci = *sci j ds psci{s)il{s)Tsci{s)e{s) . (8) 

— oo 

The dependences of fl, Tgd, and e on r and rj have been recast as dependences 
on s, as the relationship s i-^ (r, 1]) is fixed for given track parameters and 
PMT location. Figures 4 and 5 show Psci("S) for 300 MeV muons and electrons, 
obtained via Monte Carlo simulation. The distribution is relatively fiat for 
muons until the end of the track, where the ionization rate per unit track length 
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Fig. 4. Scintillation emission profile for 300 MeV muons as a function of s, the 
distance along the track. 
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Fig. 5. Scintillation emission profile for 300 MeV electrons as a function of s, the 
distance along the track. 

becomes larger. (Saturation effects are taken into account in the simulation.) 
For electrons, the distribution reflects the showering behavior of the track, 
with psci{s) rising and falling with the number of shower particles. 

We now consider an extended track emitting light with a non-trivial an- 
gular distribution, as is the case for Cherenkov radiation. We introduce 6, the 
angle to the PMT from the track, as shown in Figure 3 (right), and express 
the predicted charge /ich due to Cherenkov radiation as: 



oo 

fJ'Ch = ^Ch j dspch{s) ^{s) Tch{s) e{s) g{cos 0{s); s) 



(9) 



This expression differs from its scintillation counterpart by the presence 
of an angular emission profile g{cos6{s); s). Note that g depends on s in two 
distinct ways: the angular profile of the light changes as the track propagates 
and loses energy, and the angle 9 to the PMT changes depending on which 
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part of the track we are considering. 



Figures 6 and 7 show pch.{s) and g{cos9; s) for simulated 300 MeV muons 
and electrons. As the muon propagates, loses energy, and approaches the 
Cherenkov threshold, the rate of Cherenkov radiation per unit track length 
decreases and the Cherenkov angle becomes smaller. The scattering of the 
muon, which causes deviations of the track from its original direction, has 
been included in the Monte Carlo simulation and results in the spread of the 
angular distribution about the nominal cos^^. For electrons, pch{s) follows 
the shower profile, like Psci('S). The presence of the shower particles is readily 
apparent in the g{cos9;s) distribution, which becomes substantially wider as 
the "track" propagates. 

In writing Eq. (9), a few simplifications are made. A more rigorous treat- 
ment would involve ^ dcosOdf^ ' differential probabihty of sending an emit- 
ted Cherenkov photon in the direction {cos9,(f)), which would be integrated 
over the solid angle subtended by the PMT. By assuming that d'cosed4!'^ 
stant over the PMT face, the integration can be effected by simply multiplying 
the differential probability by the solid angle Q. Also, the azimuthal symme- 
try of the emission reduces the two-dimensional PDF to a one-dimensional 
expression, namely g'(cos^). We take g{cos9) to satisfy 



for all values of s, with the result that a factor of (An) ^ is absorbed into the 
definition $ch- 

4-3 Indirect Light (scattering, fluorescence, reflections) 

The above formalism determines the predicted charge for light arriving at 
the PMTs directly from the track without any redirection. However, the de- 
tector has sources of indirect light from scattering, fluorescence, and reflection 
as discussed in Section 3. 

The geometry for indirect light given scintillation (isotropic) emission is 
shown in Figure 8(left), where an inflnitesimal element ds along the track is 
situated at radius R from the center of the detector and at angle © relative 
to the position of the PMT. The direct light from this track element is simply 
the integrand of Eq. (8): 



1 




(10) 



-1 



d/X™ = ds $,ei Psci(s) n{s) T,ei(s) e{s) . 



(11) 
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Fig. 6. Top: Cherenkov emission profile for 300 MeV muons as a function of s, the 
distance along the track. Bottom: Two-dimensional emission profile showing the 
angular distribution of Cherenkov light as a function of s for a 300 MeV muon. 
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Fig. 7. Top: Cherenkov emission profile for 300 MeV electrons as a function of s, 
the distance along the track. Bottom: Two-dimensional emission profile showing the 
angular distribution of Cherenkov light as a function of s for a 300 MeV electron. 
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Fig. 8. Geometry of indirect light from an isotropic (left) and directional (right) 
source of light. 

An analytic expression for the indirect light would involve an elaborate inte- 
gral over emission angles and scattering points throughout the tank. Rather 
than attempt this, we observe that the value of such an integral must be 
proportional to the source strength and must otherwise depend only on the 
topological variables R and 0. The source strength dependence can be elimi- 
nated by forming a ratio of the indirect and direct light predictions: 



^sci(-R,cos6) 



indirect 
sci 



direct 
sci 



(12) 



Asci, which we refer to as the scattering table and which we build via the 
detector simulation, is a property only of the detector optics and the {R, cos 9) 
of the track element. With this table, the indirect light contribution can be 
immediately incorporated into the expression for predicted charge: 

oo 

A^sci = *sci j o?spsci(s)^(s)rsci(s)e(s)[l-F Asci (i?(s), COS e(s))] (13) 

— oo 



where the dependence of R and on s has been made explicit. 

For Cherenkov light, the situation is more complex since the light emission 
is anisotropic. Two additional variables are needed to specify the direction of 
a vector source relative to the PMT position and the tank center. We make 
the following non- unique choice for these two variables, as shown in Figure 
8(right): 

• 9: the angle between (1) the source direction vector and (2) the source-to- 
PMT ray (the same 6 as defined elsewhere). 

• (p: The angle between (1) the plane containing the tank center, the PMT, 
and the source, and (2) the plane containing the track and the tank center. 

The intensity of the indirect Cherenkov light is normalized with respect 
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to a fictitious isotropic Cherenkov source with predicted charge d/jiQ^'^^'^^" , so 
that: 

J. .indirect 

^Ch(i?,cose,cos^,0)^-/||^ (14) 



The total Cherenkov contribution to the mean predicted charge, including 
both direct and indirect components is 

.direct , .indirect 

where IjLq{^'^^ is given by Eq. (9) and /Iq^"^^^ is given by 

oo 

^indirect ^ f J^^^^^^ 

-L (16) 
X Ach (-R(s), cos6(s), cos6'(s), 4>{s))] 



Implicit in this expression is the assumption that the Cherenkov angular emis- 
sion profile does not change as the track propagates (and loses energy). This 
simplification has negligible impact, as the processes that result in indirect 
light tend to destroy the initial emission pattern anyway. 

We now have a complete charge prediction, where the total predicted 
charge for the PMT is given by: 

A* = A«Ch + fJ'sci (17) 



4-4 Computation of Integrals 



The predicted charge integrals above would be impractical to evaluate with 
sufficient spatial granularity within a maximum likelihood fit, where they must 
be evaluated multiple times for every PMT in every event. To mitigate this, 
all integrations are performed and tabulated beforehand as follows. 

The integrand of Eq. (8), which pertains to direct scintillation light, is the 
product of a "source" factor $sciPsci('S) and an "acceptance" factor 

J{s) = Q(s)T,,i(s)e(s) . (18) 



As in Eq. (8), the dependence on r and rj is recast through s. If J(s) varies 
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gradually enough, it can be approximated with a parabolic form: 

J (s) = io + iis + 32 ■ 



(19) 



An example of this is shown in Figure 9. The predicted charge then becomes 



direct 



OO CXJ CXJ 

jo j ds psci{s) + ji j ds Psci{s) s + j2 J ds Psci(s) s"^ 



(20) 



The first integral is identically unity. The remaining two depend on the energy 
Eq via Psci{s), but they depend on no other track parameters. This allows one 
to tabulate the integrals beforehand as a function of Eq, eliminating their 
evaluation from the actual likelihood calculation. 

The parabohc coefficients {ji} are obtained by evaluating Eq. (18) at three 
points along the track: the start of the track (s = 0), the midpoint of the 
Chcrenkov emission profile (s = Asmid) and twice this distance from the start 
(s = 2Asinid). The choice of where to evaluate J(s) is somewhat arbitrary; one 
could choose any three points that sample a large fraction of the range over 
which light is produced. The extension to indirect isotropic light is straight- 
forward with a redefinition of J{s) that incorporates the scattering table: 



J- 



indirect i 



s) = fl{s) Tsci(s) e(s) Asci (-R(s), cos 6(s)) 



(21) 



An analogous three-integral expression can be formed for the Cherenkov 
predicted charge: 



..direct ^ 



Ch 



jo J ds pch{s)g{cose{s);s) 

— oo 

oo 

+ ji j ds pch{s) g {cos 9 {s);s)s 

— oo 

oo 

+ j2 j ds pch{s) g {cos 9 {s);s)s^ 



(22) 



As in the scintillation case, the integrals depend on the energy of the track. 
However, they also depend on two parameters which define the PMT-track 
geometry. We choose these parameters to be the vertex-to-PMT distance, 
r(0); and the cosine of the angle-to-PMT viewed from the vertex, cos^(O). 
Labeling the integrals I^^, we obtain 

/i^r = ^Ch (ioXo°^ + ' + j22?') . (23) 
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Fig. 9. Example of J(s) for the track geometry shown in the top figure, where 
a muon 3.48 meters from a PMT moves perpendicular to the PMT axis. J(s) is 
evaluated at three points located along the track, corresponding to the three points 
in the figure. In the bottom panel, the parabolic approximation for J(s) based on 
the three evaluated points (dashed red) is compared to the exact form (solid black). 



Although the integrals {-^i"^} depend on three parameters, it is still feasible 
to tabulate their values ahead of time, eliminating their explicit evaluation 
within the likelihood calculation. 



For indirect Cherenkov light, recall that the scattering table is nor- 
malized to a fictitious isotropic Cherenkov source. Thus, the parabolic method 
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used for indirect scintillation can be applied here to give 

jindirect^^) ^ n{s) Tc^is) e{s) Ach {R{s) , COS Q (s)) , (24) 

with two energy-dependent integrals /*2\ 2 — Idspch{s)s'^ to be tabulated. 
4-5 The Time Likelihood 



The time portion of the likelihood involves the PDF /((i;x) for the mea- 
sured PMT hit time t given track parameters x. Some of the dependence on 
the track parameters can be eliminated by using the "corrected time" f^: 

f = t-to- ^(^^■^id(^o)) _ Ag^id(£;o) ^25) 



The expression removes from t three terms, namely the starting time of the 
track to, the expected time for light to propagate from the track midpoint 
to the PMT (with the speed of light in mineral oil given by c„), and the 
time for the particle to propagate from its starting point to the midpoint. 
The dependence of Asmid on the track energy Eq is made exphcit. Defining 
with respect to the track midpoint (as opposed to, say, the start of the track) 
improves the validity of the simplifications we make below. 

Given the t ^-^ substitution, we seek the PDF /jc(t?;x) for the corrected 
time given arbitrary PMT-track configurations. The space of configurations is 
five dimensional. However, producing tables of ftc {f') as a function of five pa- 
rameters is impractical. To reduce the task, we make the assumption that the 
corrected time PDF depends only on the track energy, the predicted "prompt" 
charge, and the predicted "late" charge, where "prompt" corresponds roughly 
to light arriving at the PMT directly from the track without delays induced by 
emission lifetimes, scattering, etc. Loosely, the assumption is that the shape of 
the corrected time spectrum is dominated by the physical extent of the track 
(characterized by its energy), and by the amount of prompt and late light 
reaching the PMT. The extent of the track affects the spread of possible hit 
times, since there is a spread of photon production times, while the amounts 
of prompt and late light affect the proportion of peaked "prompt" response 
to the tail of "late" response. 

This assumption reduces the configuration to three dimensions. We make 
one further simplification by assuming that for a given energy £^0; ff^if^) can 
be modeled by separate primitive prompt and late distributions, indexed by 
/^prompt and Hiate, respectively, from which we can construct the full PDF. As 
a result, the PDF is indexed not by the fundamentally two-dimensional space 
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Fig. 10. Distributions of for direct Chcrcnkov light for R5912 PMTs. The primary 
difference between the two PMT types is that the R1480 PMTs show wider prompt 
timing peaks than do the newer R5912 PMTs. 

of (/Upi.ompt,/Uiatc), but by ^Uprompt and yUlatc Separately, with the full PDF cal- 
culated on-the-fly as described below. In practice, Cherenkov and scintillation 
primitive distributions are created and used as proxies for the desired prompt 
and late distributions, respectively. 

The Cherenkov (prompt) primitive distributions are created by simulating 
particles throughout the detector with isotropically chosen directions and fixed 
energy Eq. The particles are created with direct Cherenkov light only; all other 
sources of fight (scintillation, scattering, etc.) are turned off. For each event, 
the true track parameters are used to evaluate the direct Cherenkov predicted 
charge fJ^ci^^^ for each hit. Histograms of the corrected time of the hits are 
produced for various ranges of predicted charge. Figure 10 shows three such 
distributions for 300 MeV muons. Since only direct Cherenkov ligfit is present, 
the histograms show no late-time features. The shapes of the time spectra de- 
pend on the /ich^^* values involved, with the distributions becoming narrower 
and earlier with increasing predicted charge due to the increasing probabil- 
ity that an early photon will be recorded at the PMT. The distribution in 
each predicted charge range is fit to a Gaussian distribution, and the resulting 
Gaussian parameters (mean and width) are subsequently parametrized across 
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Fig. 11. For 300 MeV muons, the mean and width parameters from Gaussian fits like 
those in Figure 10 versus Atch^^*- '^^^ dependence is parameterized as a sixth-order 
polynomial in log(/x^'£^^*). 

fj^cl^^^ using a sixth-order polynomial in log(/ich'^'^*). Figure 11 shows an exam- 
ple of these "second-level" fits for 300 MeV muons. The procedure is repeated 
at many values of Eo, with the two second-level fits providing seven parame- 
ters for the Gaussian mean and seven parameters for the Gaussian width at 
each energy. The energy dependences of these 14 parameters are then fit as 
a function of Eq to fourth-order polynomials in a third-level parametrization, 
as shown in Figure 12 for the first seven of the second-level parameters. In 
addition to conveniently summarizing the dependence of the time response on 
f^ch^^ and Eq, the parametrizations provide a smooth likelihood surface, as 
required by the minimization algorithm. 

This completes the prompt primitive distributions Gch(^"^; Eq, ij>c{^'^^)- They 
are calculated for a given predicted charge yUcT'^* ^ PMT for a track with 
energy Eq by evaluating the 14 third-level curves which parametrize the Eq 
dependence of the second-level functions. These 14 values then determine the 
two functions which parametrize the dependence of the mean and width of the 
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Fig. 12. Third-level fits parametrizing the energy dependence of the second-level 
parameters describing the Gaussian mean and width parameters. 

Gaussian parameters on ^q{^'^^ at Eq. These are evaluated at the particular 
value of Hci^^^ to determine the appropriate mean and width of the Gaussian 
distribution which describes the Cherenkov primitive distribution appropriate 
for the particular values of /XcT'^* Eq. The vahdity of the assumptions 
and simplifications can be tested directly by comparing the calculated primi- 
tive distributions with those actually observed in the Monte Carlo simulated 
events. Figure 13 shows such a comparison. 



For the scintillation (late) primitive distributions, events are generated 
with scintillation light only, and the histograms are created in bins of /U.gci^'^*. 
Each histogram is fit to a sum of two exponentials convolved with a Gaussian 
distribution. The exponential decay constants are fixed at ri = 5 ns and T2 = 
30 ns, leaving three free parameters: the time origin, the Gaussian resolution, 
and the relative weight of the two exponentials. At a given Eq, the three 
parameters are extracted in each histogram to obtain their dependence on 
nf^r^^j which is parametrized as a sixth-order polynomial in log(yUgj!['^'^*). The 
dependence of the three sets of seven parameters from the polynomials are 
parametrized as a function of Eq by fourth-order polynomials, in analogy with 
the Cherenkov case. We now have Gsci(^'^; Eq, /XgJ*'^*), the scintillation primitive 
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Fig. 13. Check of primitive distributions. The parametrized hkehhood distri- 
butions for Cherenkov hght are compared against the actual distributions for 
250 (top), 600 (middle) and 1500 (bottom) MeV muons as a function of predicted 
charge. (Note that the high-ii^O) low-fi hits, which are less well modeled, are rare.) 



distribution. Figure 14 shows the primitive distributions compared with actual 
distributions from the Monte Carlo simulation for muons at Eq = 250, 600 
and 1500 MeV. 
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Fig. 14. Check of scintillation primitive distributions. The parametrized likelihood 
distributions for scintillation light are compared against the actual distributions 
for 250 (top), 600 (middle) and 1500 (bottom) MeV muons as a function of predicted 
charge. 

To obtain the complete PDF ft<^{t'^) for a given PMT hit, we first divide the 
predicted charge into "prompt" and "late" components based on the sources: 



/^prompt 
/^late 



I .indirect i .direct i .indirect 
"T PCh "T Psci "T Psci 



(26) 
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These definitions of /^prompt and ^latc encapsulate tlie following assump- 
tions. First, all prompt light is due to direct Cherenkov light; all other light 
(including indirect Cherenkov light) follows the late scintillation-based time 
distribution. The 5% of direct Cherenkov light which is included in /iiate ac- 
counts for the late pulsing of PMTs which causes promptly arriving hght to 
appear at late times. The late time distribution is used to model indirect 
Cherenkov light since the fluorescence, scattering, and reflection processes 
which lead to indirect light have a time structure similar to that of scintil- 
lation. Further, late timing is less critical to the reconstruction, as prompt 
light provides essentially all the useful timing information. 

The /iprompt and /iiate values are used to combine the prompt and late 
primitive distributions to form the total PDF ff^it^). In the process, one must 
account for the fact that a PMT can register only one hit for a given track 
due to the latency period of the electronics. Since the time reported by the 
PMT reflects the time of the flrst photoelectron, we assume that a PMT hit 
that includes a prompt photoelectron obeys the prompt primitive distribution 
regardless of how many late photoelectrons may follow. The Poisson distri- 
bution gives us the probabilities that no prompt or late photoelectrons are 
produced based on their expected mean values: 

P(no prompt PEs) = e~''P'^°™pt ^^^^ 
P(no late PEs) = e^'^iate 



With these probabilities, one can calculate the probabihty that a given hit 
contains at least one prompt photoelectron: 

V 1 — -Pfno prompt PEs) 

Pfprompt PE present hit) = — ^ — — — (28) 

^ ^ ' ^ 1 - P(no prompt PEs)P(no late PEs) ^ ' 



We assign this as the weight Wp for the prompt primitive distribution within 
the total PDF, while wi — 1 — Wp is used as the weight for the late primitive 
distribution, yielding: 

/tc(t'^; Pq, /^prompt, /^late) = WpGchif^'-, Eq, /iprompt) + WiGsdit^] pQ, Ablate) (29) 



where Gch and Gsd are the prompt Cherenkov-based and late scintillation- 
based timing distributions, respectively. Note that Eq. (28) removes the overall 
probability of a hit occurring. Thus, even if the absolute probability of a 
prompt photoelectron is small, ~ 1 if /ii^te ^ /^prompt- 
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(a) single e/ju track 
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Fig. 15. Internal fit parameters for (a) a single muon or electron track and (b) two 
photon tracks. Each photon track includes a conversion distance parameter s. The 
two-track parameters can be constrained such that the invariant mass of the two 
photons M^^ is always M^o. 

5 The Two-track Reconstruction 



The single track model with seven parameters is sufficient for reconstruct- 
ing and z/g CCQE events as well as cosmic muons or their decay electrons. 
As discussed in Section 2, NC 7r° events require a two-track model. Figure 
15 shows the 12 parameters needed to describe two photon tracks originating 
from a common vertex. The electron track model serves double-duty as the 
photon track model, under the assumption that a showering electron is indis- 
tinguishable from a showering photon apart from the conversion distance. (The 
mean conversion length in Marcol 7 is 67 cm.) While the two-track model is 
conceptually a straightforward extension of the one-track case, complexities of 
the likelihood space require special treatment in the likelihood maximization 
algorithm. 



5.1 Two-track Charge Likelihood 

Since the charge likelihood Fq depends only on the measured charge q and 
the total predicted charge fi at each PMT, the two-track charge likelihood can 
be obtained by adding together the predicted charges from the two tracks to 
obtain the total predicted charge at each PMT: 

= /^track 1 + A'track 2 (30) 
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5.2 Two-track time likelihood 



The one-track time likelihoods account for the "first photoelectron only" 
nature of the electronics by calculating the probability that a prompt photo- 
electron is present and by then weighting the prompt primitive distribution's 
influence on the full time PDF accordingly. This scheme is extended to handle 
two tracks as follows. 

The two single-track primitive distributions, G^h ^Id formed for 
each track, where i labels the track number. In anticipation of aggregating all 
late light later, the late "scintillation" primitive distributions are averaged to 
form: 

Gsci = I (GLi + C^ci) (31) 



Of the two tracks, one will have a midpoint that is nearer to the target PMT 
than the other. This track's quantities are labeled with "n" (near); the other's, 
"f" (far). In analogy with Eq. (26), we define: 

A*prompt,n = 0.95 /^ch^n* 

/^prompt,f = 0.95 l^'^if (32) 
Ablate = A^tot A'prompt,!! /^prompt,f 

The probabilities of not obtaining a photoelectron from these sources are: 

Pn = P(no prompt PE from near track) = 6"'*?™™?*'^ 

P{ = P(no prompt PE from far track) = 6"''?'"°™?*'^ (33) 

Pi = P(no late PEs) = e'/'iate 

from which the following weights are obtained: 



1 - PnPfPi 
1-Pf 



1-PfPi ' 



Wl — 1 — Wn — Wl 
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where Wn is the probabihty that a prompt photoelectron from the near track 
exists given that any photoelectron exists, etc. The weights are used to com- 
bine the two prompt primitive distributions and the averaged late primitive 
distribution: 

Mt') = Gch,„(t^) + w, Gci^An + ^1 Gsciin (35) 



resulting in the complete two-track corrected time PDF. 



6 Minimization of — log(£) 
6.1 One-track 



With the likelihood C defined, the parameter set x that minimizes its neg- 
ative logarithm F = — log(£) must be found. In the single track minimization 
process, two issues arise. 

First, the energy parameter Eq is tied to the geometry of the event via 
the track profiles p{s;Eo). If the spatial parameters (.tq, etc) are varied si- 
multaneously with Eq, the minimization algorithm can get confused by this 
correlation. The solution is to iterate the minimization process, as described 
below. 

A second issue arises from the discrete PMT lattice, which imprints a 
small fluctuating signal on i2(x) despite the smooth parameterizations of the 
input tables. Minimization algorithms that rely on gradients or that do not 
have controllable step sizes are thus poorly suited to the problem. Therefore, 
the SIMPLEX method in Minuit[ll] is chosen over the MIGRAD method. 

To initiate the minimization process, a vector of seed parameters is ob- 
tained from a fast fitter. In the first iteration of the minimization, the energy 
Eq is held at its seeded value while the six remaining parameters (xq. Ho, Zo,to, 
9o, (f)o) are varied via Minuit /SIMPLEX to find a temporary minimum of F. 
These six parameters are then fixed, while Eq is freed for a second iteration. 
Finally, Eq is fixed once more with the six other parameters varied (as in the 
initial iteration), to find the final minimum of F. The parameters from the 
final SIMPLEX call are returned as the best-fit set x. 
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Fig. 16. The two it fitter configurations discussed in tlie text. 
6.2 Two-track 

Figure 16 shows two situations that may arise in fitting a vr'' event. The 
event has two photon tracks whose Cherenkov rings are represented in black, 
while green and red represent two ring configurations corresponding to dif- 
ferent starting values for the parameters. In case (a), the 12 parameters are 
near the correct configuration. In case (b), the parameters are such that both 
tracks are directed toward the dominant ring. The latter case corresponds to a 
local minimum which may offer a worse likelihood than case (a); however, the 
minimization algorithm may be trapped by the fact that any small changes 
to the parameters result in an increased value of F. For example, sweeping 
one track over to the smaller black ring involves passing through a region with 
little detected Cherenkov light. The intermediate states form a barrier of dis- 
favored values of F that the minimization algorithm code will have difficulty 
overcoming, especially if several parameters must be adjusted simultaneously 
even to make the distant configuration an improvement. 

The two-track fits require a minimization approach that avoids these traps 
in the likelihood surface. Monte Carlo 7r° events were used to identify both 
common and rare trapping scenarios, and the minimization algorithm ad- 
dresses each of these through the following collection of seed configurations 
and minimization sequences:n 
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• The conversion lengths Si and S2 are seeded with either 50 cm or 250 cm, 
leading to four possible pairs. 

• 6i and 0i are seeded with the results from the electron single-track fit or 
with one of the eight perturbations (see below). 

• $2 and 02 are seeded with the "best" directions from a full grid of tested 
directions (see below). 

• The four-vertex of the event is seeded with the four-vertex from the electron 
single-track fit shifted according to Si, 6i and 0i. 

• El is seeded with approximately Eq from the electron single-track fit, while 
the seed for E2 is, when possible, the energy needed to give an invariant 
mass equal to M^o. 

For the last item, the energy seeds are based on an empirical expression 
that accounts for the second ring's energy contribution to the one-track fit 
energy as a function of the angle between the two 7 tracks. The function is 
obtained by simulating NC tt" events and reconstructing them with the single 
track electron fit. Using the true photon energies {Ei > E2) and opening angle 
9.y.y, the relationship between {Eq — Ei)/E2 and cos 6'-^-^ is parametrized by a 
linear function. This function is used to determine the seed Ei and E2 values 
for a given Eq and cos6'^^. 

The nine (1+8) possible seed directions for track 1 are created as follows. 
The spatial distribution of PMT charge is projected onto a plane perpendicu- 
lar to the electron fit direction. The major axis of the covariance ellipse of the 
resulting two-dimensional distribution is found. The nine possible seed direc- 
tions are rotations of the best electron-fit direction by 0, ±0.159, ±0.450 and 
±0.644 radians parallel to the major axis and ±0.159 radians perpendicular 
to the major axis. 

The seed directions for track 2 come from either a "coarse" 24x12 {6, </>) 
grid or a "fine" 50x25 {6, (p) grid. The grid type is specified at run time. 

Of the 9x(24xl2)=2592 or 9x (50x25)=11250 possible combinations of 
track 1 and track 2 seed directions, only six are used to seed actual minimiza- 
tion sequences. They are the six combinations with (1) the best total likelihood 
(time and charge), (2) the best charge likelihood when the estimated track en- 
ergie^ are similar (0.5 < ^ < 2), (3) the best charge likelihood when the 
estimated track energies are dissimilar, and (4) — (6) the second-best likelihood 
in each of these three cases. A small exception: combinations involving one of 
the eight perturbed track 1 directions and having estimated track energies that 
differ by more than a factor of five are not considered. In these energetically 
asymmetric cases, the electron fit direction would be a good representation of 
track 1 since the influence of track 2 would be small, so there is no reason to 



Recall the energy seed procedure discussed above. 
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consider alternate track 1 directions. 



The six best and second-best track 1 and track 2 direction and energy 
combinations are chosen anew for each permutation of (si,S2), leaving 24 
complete parameter sets in total. Each set seeds two Minuit sequences. In 
sequence 1, the SIMPLEX minimization routine is first run with the energy 
and angle parameters held fixed; once a minimum is found, SIMPLEX is called 
again with all parameters free. In sequence 2, only the energy parameters 
are held fixed in the initial SIMPLEX call. This results in 48 fit sequences, 
each of which ends with a completely free SIMPLEX minimization. The best 
parameter set {i.e., the one that results in the lowest F) from the forty-eight 
final SIMPLEX calls is reported as the best parameter set x for the two-track 
fit. During all minimizations, the photon conversion points are constrained to 
be within the main detector volume. This prevents the fit from removing the 
influence of a track through an inflated conversion length. The event vertex 
itself, however, is not constrained. 

The SIMPLEX minimizations can be performed with a constraint on the 
invariant mass. This is accomplished by removing E2 as a free parameter and 
setting it via the relation: 

= 2E^{1-Lse,,) ^^^^ 



where 6^^ is the angle between the photon tracks. This fixed-mass mode is the 
actual 7r° hypothesis, while the free-mass mode allows for mass reconstruction. 
Both are used for 7r° identification, with the former lending its maximized 
likelihood Cj^o and the latter providing a reconstructed mass M^^. 



7 Performance 



The performance of the reconstruction algorithm on simulated neutrino 
interactions (i^^ CCQE, CCQE and NC vr" events) generated with NU- 
ANCE [8] and propagated through the detector Monte Carlo simulation is 
shown in Figures 17— 2l|^The electron and muon one-track fits are applied 
to the and CCQE events, respectively, while the two-track fits, with 
and without the mass constraint, are applied to the NC 7r° events. A minimal 
selection is applied to ensure that the events are contained within the main 



^ For comparisons of Monte Carlo to data and for discussion of the development 
and vahdation of the Monte Carlo, see Refs. [9] and [10]. 
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detector region and, for plots not directly related to vertex reconstruction, 
within the fiducial volume of the analysis. 

A few observations about the figures, particularly Figure 18: (1) The reso- 
lutions for z/^ CCQE events degrade as the energy approaches the Cherenkov 
threshold from above, as most of the information in the event comes from the 
(rapidly decreasing) prompt Cherenkov light. (2) A low-energy vr^ offers little 
information on its direction, as all that is observed is the isotropically produced 
77 final state. Thus, the angular resolution (but not the vertex resolution) be- 
comes poor. (3) The vertex resolution for vr'' events is worse than Vf. CCQE 
events since the former has spatially displaced light production thanks to the 
~70cm average photon conversion length. (4) Typical performance numbers 
for z/e CCQE events: 20 cm vertex resolution, 12% kinetic energy resolution, 
and 4° angular resolution. (5) Typical performance numbers for z/^ CCQE 
events: 10 cm vertex resolution, 8% kinetic energy resolution and 2° angular 
resolution. (6) Figure 19 shows a 7r° mass resolution of ~13%. 
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Fig. 17. Reconstructed radial event position (left) and kinetic energy (right) plotted 
against their true values. The top, middle, and bottom panels show Vp, CCQE, 
CCQE, and NC tt^ events. 
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Fig. 18. Radius, kinetic energy, and direction resolutions. The jitter is due to limited 
test sample statistics in some regions. The V(, curves cut off at 100 MeV due to an 
unrelated upstream cut on the Monte Carlo sample. 
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Fig. 19. Fitted invariant mass M^^ for Monte Carlo-simulated NC vr'' events in three 
low energy regions (left) and three high energy regions (right). 



NC n" NC n° 




against true values. The left (right) panel shows the higher (lower) energy 7 from 
each event. The association of each fitted track to the underlying true 7's is, in 
general, ambiguous. For these plots, we choose the assignment that gives the smaller 
combined energy and direction discrepancy. 
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Fig. 21. Energy (left) and direction (right) resolutions for 7's in simulated NC vr^ 
events. For a given true 7 energy, a dominant 7 has expectedly better resolution 
than does a subordinate 7. 

8 Particle Identification 



The maximum likelihoods returned from each fit can be used for hypothesis 
testing. The two quantities 

7^e/^ = log = log £e - log (37) 



and 



7^e/7r = log = log Ce - log £^ 



(38) 



are used to determine for a given event whether the electron hypothesis is 
preferred over the muon and 7r° hypotheses. In these expressions, Ce, C^ and 
are the maximized likelihoods returned by the electron, muon, and (fixed- 
mass) two-track fits, respectively. 



8.1 Electron/ Muon separation 



Figure 22 shows for simulated Ve CCQE events and CCQE events the 
distribution of T^e/^t as a function of the energy reconstructed by the electron 
fit. For each sample, the events have been subject to a set of preselection 
criteria used in the Ve appearance analysis. These require that there is only one 
time-cluster of PMT hits in the event, eliminating obvious CC events that 
produce decay electrons. Cosmic backgrounds are eliminated by requiring the 
minimum number of PMT hits in the main region to be >200 and the number 
of veto hits to be <6. Furthermore, the average time of hits is required to lie 
within the expected beam delivery window from the Booster. 
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Fig. 22. Distribution of IZ^j ^=\og{Ce/ C^) for Monte Carlo simulated Ve CCQE 
events (top) and CCQE events (bottom) as a function of reconstructed energy 
(from the electron hypothesis fit). 
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Fig. 23. Left: Distribution of 7^e/^= log(£e/>C7r) for Monte Carlo simulated Ve CCQE 
events (top) and NC vr'' events (bottom). Right: same for M^^. The black lines on 
each distribution show the selection used in the MiniBooNE i/^ — > i^e oscillation 
analysis [4]. 

For the z/g CCQE events, 7le/f_i tends to take positive values, indicating 
that the fit to the event with the electron hypothesis is favored over the muon 
hypothesis. Likewise, TZe/fj. tends to be negative for the i/^ CCQE events, in- 
dicating that the muon hypothesis fits these events better than the electron 
hypothesis. At high energies (>1 GeV), the electron/muon separation is aided 
by the fact that the muon pathlength grows approximately linearly with en- 
ergy, while the electron shower grows more slowly. 



8.2 Electron/iT^ separation 



The left plots in Figure 23 show the distributions of TZ function of 

energy reconstructed by the single-track electron reconstruction for simulated 
Ue CCQE events and NC 7r° events. The quantity TZe/n tests whether a given 
event fits better as a single electron track or as a 7r°. Overall, electron/vr^ 
separation becomes more difficult at high energies as the energies of the decay 
photons in 7r° events become more asymmetric and the shower fluctuations 
in a single electron event get large enough to mimic the presence of a second 
photon. Also, 7r° events in which one of the two photons leaves the detector 
unconverted present an essentially irreducible background at low energy. 
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Another quantity that can be used for electron/vr" separation is My-y, the 
invariant mass of the two photons returned by the free-mass two-track fit. The 
NC tt" events peak at the vr*^ mass, whereas the z/g CCQE events peak toward 
zero. This is seen in the right plots in Figure 23, where the M^^ distribution is 
shown as a function of reconstructed energy. The use of the true 7r° mass in the 
fit's seeding procedure (Section 6.2) induces no appreciable high-mass peak in 
the z/g CCQE events while helping considerably in the correct identification of 
NC 7r° events. 

The Pe appearance analysis at MiniBooNE utilizes both M^^ and TZe/-K to 
separate z/g CCQE events from NC vr''. 



8.3 Particle Identification Performance 



The background rejection levels reachable with these fitter-derived quan- 
tities are shown for Monte Carlo events in Figures 24 and 25. Figure 24 shows 
the efficiency, after fiducial volume and cosmic rejection cuts, for selecting 
signal Ve CCQE events alongside the efficiency (or misidentification rate) for 
Vfj, CCQE events. Three prototype selections are used, each chosen to give an 
approximately fixed signal efficiency at all energies. Note that the misidentifi- 
cation rate is shown with a logarithmic vertical scale. Figure 25 shows similar 
plots for 7r° rejection, one for the TZ^/t^ cut and one for the M^^ cut. 
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Fig. 24. Performance of three prototype TZf.i ^=\og{Ce/ Cfj) cuts for the z^e selection. 
Left: Signal (fg CCQE) efficiency versus the electron reconstruction's energy. Right: 
Misidentification rate of CCQE events versus the electron reconstruction's energy. 

Misidentification rates of a few percent for z/^ CCQE events are seen for sig- 
nal efficiencies from 70 — 90%. (Note that z/^ CCQE events are further reduced 
by about an order of magnitude by recognizing the delayed muon decay.) The 
two 7r° variables indicate that the misidentification rate increases with energy; 
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Fig. 25. Performance of three prototype 7^6/.^-= log(£e/^7r) (top) and M^^ (bot- 
tom) cuts for the Ue selection. Left: Signal {ve CCQE) efficiency versus the electron 
reconstruction's energy. Right: Misidentification rate of NC vr'^ events versus the 
electron reconstruction's energy. 



this is expected, since the two decay photons become highly asymmetric ener- 
getically and since the Cherenkov rings can achieve a greater degree of overlap. 
The rise in the 7r° misidentification rate at low energies is due to events with 
an escaping photon (corresponding to the cluster of ti^ events with positive 
7^e/^ at 50-200 MeV in Figure 23). Between 200 MeV and 600 MeV, where 
most of the NC tt*^ events lie, a misidentification rate of less than 1 percent can 
be achieved with 40% efficiency using only the 7^e/7r cut, while rates of a few 
percent can be achieved with only the M^^ cuts. In the MiniBooNE — > v^, 
oscillation analysis, 50% signal efficiency and 1% vr" misidentification rate are 
obtained by requiring both M^^ < M^^''{E) and TZe/-n > Tl'^/'^{E), where the 
"max" functions are quadratic in the energy E obtained from the single-track 
electron fit. These quadratic selection boundaries are indicated by the solid 
black lines in Figure 23. These cuts, along with an analogous TZe/fi cut, were 
tuned to optimize sensitivity to LSND-like oscillations. 
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9 Summary 



A maximum likelihood reconstruction algorithm is used at MiniBooNE to 
determine track parameters under electron, muon, and 7r° hypotheses. Under- 
lying the likelihood is a track and detector model that calculates for a given 
parameter set the charge and time PDFs expected for each PMT, accounting 
for the spatially extended production of Cherenkov and scintillation light as 
well as the effects of indirect light from subsequent optical processes. 

The maximized likelihoods for a given event obtained under different event 
hypotheses are used for event selection. In particular, the ratio of the likeli- 
hoods under the electron and muon models is used to suppress z/^ CC events in 
the MiniBooNE oscillation search; the ratio of the likelihoods under 

the electron and two-track vr*' fit (with fixed invariant mass) and the invariant 
mass obtained from the free-mass two-track fit are used to suppress NC 7r° 
events. While the reconstruction has been developed within the context of the 
MiniBooNE — > Ve oscillation search, other experiments employing similar 
Cherenkov detection techniques should be able to use the techniques discussed 
here. 
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